Bloch Walls and Macroscopic String States in Bethe's solution of the Heisenberg 

Ferromagnetic Linear Chain 
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We present a calculation of the lowest excited states of the Heisenberg ferromagnet in 1-d for any wave vector. These 
turn out to be string solutions of Bethe's equations with a macroscopic number of particles in them. These are identified as 
generalized quantum Bloch wall states, and a simple physical picture provided for the same. 
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Introduction The question of elementariness of excita- 
tions in low dimensional magnetic systems is receiving 
much attention currently. It was perhaps first addressed 
in the seminal work of Bcthe in 1931 (T). In addition to 
providing the celebrated Ansatz named after him, Bethe 
asked if Bloch's magnons are the "most elementary" ex- 
citations in 1-d. He came to the conclusion that they 
were not, and instead found that the bound states of spin 
reversals were. After the original paper of Bethe, the fer- 
romagnet has received comparatively less attention 
than its antipode, namely the antiferromagnet One 
source of revival of interest in the ferromagnet is in con- 
nection with stochastic dynamical systems, albeit with 
a complex Aharonov Bohm magnetic flux [QD ■ Another 
notable recent exception is a work by Sutherland Q , who 
shows that the excited states of the ferromagnet contain 
a singlet state at momentum tt, with an excitation energy 
(EE) that is very low, of 0(1/N), where N is the length 
of the ring. 

At the semiclassical level, domain wall arguments 
lead one to expect in dimensions d ( with volume = N ) 
the Bloch wall excitations to be of 0(N d ~ 2 ) , and hence 
to be amongst contenders for the lowest EE in d = 1. 
Such "large deviation" excitations carry spin as well as 
momentum as we show below. These configurations of 
spins will be discussed within the context of Bethe's 
Ansatz (BA) for the s = 1/2 ferromagnet presently. 

In this work we ask ( and answer) the following ques- 
tion: For a given value of the total momentum or total 
spin of the Bethe ferromagnet, what is the lowest excited 
state? The nontriviality of the question arises from the 
fact that within the famous Bethe formula for the bound- 
state of n magnons, <jJBethe(q) — — cos (9))i the lower 
limit on the total momentum q depends implicitly upon 
n. Its dependence has not been fully explicated earlier, at 
least as far as we could find in the literature. In this work, 
we use a combination of exact diagonalization and ana- 
lytic methods to attack the problem. A new and essential 
tool that we develop is the argument of continuity of cer- 
tain regular root solutions with respect to the density, 



regarded as a continuous variable, leading to a differen- 
tial equation formulation of Bethe's equations (BE) that 
bypasses the knowledge ( or otherwise ) of the quantum 
numbers. 

Our findings are readily stated: The lowest excitations 
for any momentum q arise from special string solutions 
of BE. These special string solutions involve a macro- 
scopic number of particles in a given string, and hence 
we call them macroscopic strings. Sutherland's solution 
at momentum tt is a particular case of these. We have 
found the lowest solution for every value of the total 
momentum, these correspond to a definite value of to- 
tal spin as well. Such states are of the type that one 
would expect from Bethe's formula for n magnon bound 
states, with n of the size of the lattice. The formula 
of Bethe cannot, however, be used for such large bound 
states since we show that there are significant corrections 
to the traditional assumption of a uniformly spaced ver- 
tical Bethe string in the complex plane: the curvature 
and nonuniformity of spacing produces essential differ- 
ences. Our states can be represented by a new formula: 
uj BW (q) = ^Jq(l - q/{2ii)) for < q < tt. We find 
that the solution for a given q corresponds to a particu- 
lar value of total spin, Stot = N(l/2 — q/2n). If we write 
n the spin deviation from the saturated ferromagnet ( 
S* = N/2 — n) in terms of the density d as n = d N, the 
spectrum can be written as lobw — ~^~Jd(l — d) with 
q = 2nd and Stot = N(l/2 — d). We show finally that 
these states correspond to generalized quantum Bloch 
walls [Jul: thus the states with lowest EE at any wave 
vector of 0{\) are Bloch walls fli}] . 

Before presenting our calculations, we also note that 
the bound states of Bethe for s — 1/2 have been iden- 
tified recently with solitonic excitations || of the non 
linear classical, i.e. large s, Landau Lifschitz equations 

S{x) — S(x) x d 2 / dx 2 S(x). Several explicit solutions of 
these are known pfl| , pT| |. We note that Bloch walls form 
a certain class of exact solutions [|]Jll]]of the nonlinear 
classical equations , namely non linear spin waves. 
Calculations: We consider the ferromagnetic Heisen- 
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berg Hamiltonian: H = 2J^ Z=1 N [s 2 — St ■ Si+%], where 

for most part we consider s = 1/2 and Si = i<7j, in terms 
of the usual Pauli spin operators. 

Our first result is from observations on wavefunctions 
of the energy eigenstates obtained from exact numerical 
diagonalization of chains upto length N = 16 in a mo- 
mentum resolved basis. We observed that in all cases, 
the lowest state at momentum q = 2nn/N has total an- 
gular momentum S = N/2 — n. This implies that this 
state can be obtained in the n particle sector where it is a 
"maximal" , i.e. highest weight state. The corresponding 
Bethe wavefunction has all pseudo momenta nonvanish- 
ing, and the momentum and particle density related as 
q = 2nd. 

We next write the BE for the Heisenberg chain in the 
Orbach parametrization as: 



Nf(ai) = 2nh + £ /(( 



0/2) 



I = 1,2. ..71 (1) 



where f{x) = ±log(f±f) = 2ArcCot(x), a t = cot(*i/2) 
and {/;} are the Bethe integers, and ki the Bethe pseudo 
momenta. We take the branch cut of f(x) to be on 
the imaginary axis running from — i to +i. The en- 
ergy and total momentum are given respectively by: 

e = /Eti,ntHj and 1 = W Ei=i,n^i- Let us note 
that we are interested in the lowest energy states for a 
given q, requiring a knowledge of the integers /; . These 
integers, as shown by Bethe, differ by 2 for scattering 
states, and by either one or zero for bound states in gen- 
eral. Eliminating the integer sets with zeros in them, a 
very plausible state is one with Ii = 1 for 1 < I < n, and 
indeed we found from numerical studies of BE for small 
N and small n ( with n << N), that this was indeed 
so: the resulting state is invariably the lowest energy 
state for small q. Emboldened by this exercise we found 
the following exact solution analytically in the limit of a 
thermodynamic n as well as N, but at low density i.e. 
d = n/N << 1 . Since we are interested in excitation en- 
ergies that are vanishing in the thermodynamic limit, the 
corresponding variables ai scale with system size and it is 
convenient to introduce new scaled variables zi = oti/n. 
In terms of these, the lhs of Eq(|l|) becomes 2/{zid) on us- 
ing the large x expansion of f(x) and ignoring terms of 
size 0(1/N 2 ). On the rhs we cannot make the expansion 
in general since, at high densities, a core is formed ||] 
where the separations n(zt — z m )/2 between pairs of par- 
ticles can become arbitrarily close to the branch points 
of f(x). In the low density limit we find it is possible to 
obtain a perturbative solution, using the crucial obser- 
vation that the typical interparticle separation ~ l/Vd, 
and hence we can use the smallness of d as an expansion 
parameter in a perturbative sense. Setting // = 1 we find 
the approximate equation: 



1 , 2d v-^ 
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Note that Eq(||) has corrections from the expansion of 
the phase shift, that is typically of 0(d) smaller than the 
least term retained. The solution of the system Eq(||) 
can actually be found exactly jlj|, but we save it for a 
future publication. We find at low densitites the following 
result: 
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(x 2 + 1 - -) + 0(Vd), (3) 



where xi satisfy H n (^/nxi) — 0, H n being the nth order 
Hermite polynomial. In the limit of large n the Xj form a 
continuum stretching from — \[2 to \[2 with the familiar 
semicircular density of states p(x) = — v2 — x 2 . This 
solution can be used to obtain the energy to order d 2 . 
The energy Ne — 4J/(nd) J2i 1 1 z 1 for low d can be found 
as Ne = 4Jn 2 [d + d 2 {3 < (/?f } ) 2 > -2 < 0<{ 2) >} + 
0(d 3 )], where the averages are normalized sums over the 
indicated variables. Using the explicit expression Eq(|^) 
and converting the sums to integrals over the semicircular 
density of states we get finally the low density formula: 
Ne = AJn 2 d(\ - d) + 0(d 3 ). Below we will argue that 
there are no corrections to the above formula beyond the 
first term: it is exact!. Thus provisionally we write 



<l 



Ne = 4:Jn 2 d(l - d) = 2Jnq(l - — ) 

2tt 



(4) 



We note that the low density Equation(||) must be aban- 
doned once the minimum separation n{Z{ — Zj) hits the 

value 2i, this happens at 2 = \n(zo — zi)\ — ~r^[~ n \{ x o ~ 
xi). However n\(x — £Ci)| = l/p(0) = 7r/v2, thus 
d ~ 1/4. Indeed we found for small systems that d > 1/4 
cannot be treated easily numerically: the new difficulty 
is that the quantum numbers are no longer simple as we 
discuss below. For d < 1/4 the low density result Eq(|^) 
and the full solution of Eq(Q) |l5| are extremely close. 

We now discuss the techniques used for solving Eq(Q) 
numerically at larger densities. There are two main prob- 
lems. One is that the integers jump around in a compli- 
cated way that is not known beforehand and it is clearly 
not feasible to try all combinations. The second problem 
is the formation of the core i.e. successive roots that 
are placed very close to a separation 2i, which causes sin- 
gularities in the equation and results in numerical inaccu- 
racies. Our strategy is to start from low densities where 
we know the roots, and change d slowly and study the 
evolution of the roots. For this we first convert Eq(|l|) into 
a set of first order ordinary differential equations (ODE) 
with \og(d) as a time like flow-parameter. Taking the 
derivative of Eq([j]) with respect to d and defining new 
variables t = logd and /; = z;e* we obtain: 
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where A; m = 



(5) 



2e*/i 



t^2 



Ti [4e 2t /n 2 + (/, - / m ) 2 
1 - ^ A Zm 



m ^ I 



and the derivatives are with respect to the "time" vari- 
able t. The time derivatives of the integers are delta 
functions and hence drop out of the equations at almost 
all times. Also since the flow of the roots themselves is 
smooth, we expect the delta-function singularities to be 
precisely cancelled by other terms in the equation. Hence 
in evolving the above ODE we can drop the first term 
on the right hand side of the equation. This immedi- 
ately solves the problem of our lack of knowledge of the 
integers since they do not occur anywhere else in the dif- 
ferential equations. From these solutions we can recover 
the integers and use them in the root-finder to get more 
accurate solutions. We find that till densities around 
d = 0.45 the solutions obtained from the ODE are very 
accurate. At higher densities, significant numerical er- 
rors show up because of the singularities associated with 
the core. In these cases we correct our ODE solutions 
by fixing the core by hand and using the root-finder to 
self-consistently solve for the roots outside the core. We 
plot in Fig. |l| the solutions, for a system of 16 particles, 
at different densities. The table below shows the integer 
sets at four different densities (since they occur in pairs 
we show only half of them) : 
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We now discuss the energies that we obtain from the 
BA solutions. For N < 16, we have verified that all the 
solutions obtained from the numerical solutions of the 
BE using the above scheme, match with those obtained 
from exact numerical diagonalization. With the BE we 
can go to much larger system sizes. We find that the gap 
vanishes at every finite q, with system size dependence 
~ 1/N. In Fig. H we plot the system size dependence 
of the gap at two densities, namely at quarter and half 
fillings. The latter case corresponds to the q — tt state 
considered by Sutherland and we verify his result NSE = 
Jir 2 . At d — 1/4 NSE seems to asymptote to the value 
3J7r 2 /4. Remarkably both these asymptotic values of 
NSE at densities 1/2 and 1/4 can be obtained from the 
low density formula in Eq([|) . In Fig. || we also plot the 
energy-wavevector curve, obtained from the solution of 
the 16-particle problem and compare it with Eq(^). Note 
that the discrepancies at large q are finite size effects and 
would vanish in the ./V — > oo limit. 



Variational results Having found the excited states, 
we now turn to the explicit connection with Bloch wall 
states. We now show, remarkably enough, that the ex- 
pression Eq(|) for the gap can be obtained from a simple 
variational calculation. We work with arbitrary spin s of 
the particles. A neat way to generate Bloch walls is via a 
unitary rotation operator [jljj Q — exp i^- J2 m m ^m ac ^- 
ing upon an appropriate state, |0„ >= (S' ( ^) n |/erro >, 
where S~ = J2j exp(iq.rj)Sj~ is a spin wave creation 
operator carrying momentum q and \ferro > is the 
state with all spins up. Using [H, Sq] — 0, we see 
that |0„ > is the ground state in the n-particle sector 
with zero total momentum. Thus finally we write the 
variational Bloch Wall state \BW >= Q\0 n >■ Us- 
ing the quasi commutator QS~ — S~ +27r , N Q, we find 

that \BW >— (S 2ii ^ N ) n \ferro >, and has total momen- 
tum q = 2nd. In the semiclassical limit, s » 1, the 
above state is readily visualized as classical spins that 
are tipped from the z axis and rotate along a cone. The 
variational calculation of the excitation energy SE, i.e. 
< O n |QT7^Q|o n > can be done easily by transforming the 
rotation onto the spin operators. Writing n = dN, with 
< d < 2s, we find SE = ^-d(2s ~ d). At s = 1/2 this 
is also the result of the calculation in Eq(Q). 

At this point we admit that we were surprised, as 
the reader might well be, that the results of an elabo- 
rate bound state calculation with a macroscopic num- 
ber of complex roots agrees with the result of a sim- 
ple looking variational wavefunction that resembles a 
Bose condensate of spin waves. This phenomenon is 
presumably a consequence of the shallow nature of the 
bound state. We note that the variational states satisfy 
V< ft 2 > - <H > 2 / < H 0(1/ VN) which shows 
that in the N — ► oo limit these become exact eigenstates, 
thereby providing independent evidence for the exactness 
of the main result of our work. 

For large s, the energy as well as momentum of these 
states agrees with the semiclassical estimates |ll[ using a 
Poisson bracket structure to construct a semiclassical mo- 
mentum operator. The variational results for all values 
of s thus collapse onto the same formula. The conjecture 
of ||To| , [rif implies just this kind of a result, but for the 
solitons, i.e. for bound states with small number of spin 
deviations. Needless to say, we believe that our varia- 
tional results are exact (in the thermodynamic limit) for 
all spin, since we have established them at s = 1/2 and 
also for very large s. 

We note that the Bloch wall states \BW > carry a total 
spin that is easy to calculate using simple extension of 
the above calculation: < Sf ot >=< S z > 2 = N 2 (s - d) 2 . 
Finally we note that these variational results for Bloch 
walls are generalizable to higher dimensions. 

Concluding remarks We finally note in summary that 
the Bethe formula for n magnon bound states: uJBethe — 
J 2.(1 — cos(g)) is (a) exact for q >> 2irn/N, (b) invalid 
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for q < 2im/N and (c) has significant corrections when 
q r~j 2irn/N The evidence for (a) is in Bethe's paper 
itself, the corrections to it are exponential in N, (b) is 
numerical. Evidence for (c) has been presented in this 
paper, where we find instead: ujbw = 2irJ /N\q\(l — J^). 

1 Also at the Poornaprajna Institute, Bangalore. 
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been used earlier. It is crucial in circumventing the prob- 
lem of lack of knowledge of the Bethe integers. 

[17] The Q operator is used extensively in studies of the an- 
tiferromagnet following E Lieb, T Schulz and D Mattis, 
Ann. Phys. 16, 407 ( 1961). For ferromagnets see 

[18] Note that this is true even for small n. For example for 
n — 2, q — 4tt/N the roots, from Eq(^), are at ai.2 
Y - ± i\/~N /ty instead of Bethe's N/tt ± i. As noted by 
Bethe, for the 2-magnon case the string solution holds 
for q ~ 1/^/N. 
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FIG. 1. The figure shows the Bethe curves for a 16-particle 
system at different densities (values indicated along each 
curve). The inset compares the numerically obtained roots 
(points), for a 32-particle system at d = 0.2, with the exact 
low-density result (solid line). 
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FIG. 2. The gaps at different wave vectors as obtained 
from the solution of the BA equations for 16 particles. The 
gap as given by Eq(Q) is also shown. The inset shows depen- 
dence of the gap on system size at two different wave vectors. 
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